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Preface 


This  report  represents  the  combined  results  of  the  authors  individual  investigations.  Robert 
Mailloux  conceived  the  project  and  completed  the  analysis  of  Methods  1  and  2  as  well  as  formulas  for 
all  the  basic  array  parameters.  Edward  Cohen  performed  the  analysis  of  Method  3.  wrote  and 
debugged  all  computer  programs  and  derived  the  asymptotic  dependence  in  the  Appendix. 
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Statistically  Thinned  Arrays  with 
Quantized  Amplitude  Weights 


1.  INTRODUCTION 

Skolnik,  et  al1.  and  others2,3  have  Investigated  the  use  of  statistically  designed  density  tapers  to 
control  antenna  pattern  sidelobes.  In  the  procedure  described  by  Skolnik,  all  of  the  elements  of  the 
array  were  excited  with  the  same  weights  and  the  density  of  elements  was  made  proportional  to  the 
amplitude  of  the  aperture  illumination  of  a  conventional  "niled”  array.  The  resulting  ensemble- 
averaged  power  pattern  was  shown  to  be  equal  to  the  sum  of  the  filled  array  pattern  and  a  distribution 
that  was  not  angle-dependent,  and  which  represented  an  average  sidelobe  level.  A  recent  paper  by 
Numazaki,  et  al^  combines  statistical  thinning  with  a  deterministic  procedure  to  achieve  improved 
patterns  in  pre-selected  angular  planes  (the  principal  orthogonal  planes  and  the  diagonal  planes). 
The  authors  also  extended  their  procedure  to  address  multiple  quantized  weights.  The  present  work 
also  addresses  thinned  arrays  with  multiple  quantized  weights,  but  is  primarily  concerned  with 
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evaluating  the  ensemble- average  array  parameters  for  several  quantizing  algorithms.  This  study  was 
motivated  by  the  fact  that  if  array  amplitude  weights  are  applied  by  solid  state  modules,  then  it  is 
often  practical  and  more  economical  to  select  modules  with  a  small  number  of  discrete,  or  quantized, 
amplitude  levels,  and  to  organize  the  array  into  annular  rings  with  the  quantized  levels.  Since  an 
array  with  quantized  amplitudes  may  produce  high  sidelobes,  the  goal  of  this  study  was  to  use 
statistical  thinning  to  smooth  the  average  amplitude  illumination  and  reduce  peak  sidelobe  levels. 
The  changes  in  average  sidelobe  level  gain  and  effective  isotropic  radiated  power  due  to  this  design 
procedure  were  then  evaluated. 

2.  METHODS  OF  STATISTICAL  THINNING  AND  QUANTIZING 

Several  design  procedures  will  be  discussed  that  require  an  "ideal"  continuous  distribution  as  a 
starting  point  to  synthesize  the  desired  arrays.  Figure  la  shows  a  representation  of  a  continuous 
amplitude  distribution  that  will  be  considered  an  ideal  illumination  designed  to  produce  a  low 
sidelobe  pattern  circular  aperture.  Throughout  this  report  the  aperture  distribution  will  be 
represented  by  the  circular  Taylor  patterns.5  The  circular  aperture  illumination  will  be  sampled  by  a 
rectangular  grid  array  (Figure  lc)  with  elements  located  0.5X  apart,  as  done  in  the  earlier  Skolnik 
study.  Figure  lb  shows  a  step-wise  continuous  aperture  distribution  composed  of  three  levels  that 
define  the  set  of  available  weights.  By  proper  selection  of  the  levels,  which  we  refer  to  as  amplitude  or 
weight  quantizations,  and  applying  one  of  the  three  thinning  procedures  described  below,  one  can 
approximate  the  ideal  aperture  illumination  (A^  and  achieve  low  sidelobe  average  radiation  patterns. 
Throughout  much  of  this  analysis  we  shall  use  Skolnik's  notation,  and  denote  aperture  sums  as  one¬ 
dimensional  sums  over  the  index  n  while  understanding  this  as  a  shorthand  for  the  conventional  two 
dimensional  array  summations. 

The  three  thinning  methods  used  in  this  study  are  depicted  in  Figure  2.  The  ideal  aperture  taper, 
indicated  dashed,  is  approximated  either  by  omitting  elements  entirely  (thinning)  or  by  selecting 
weight  values  from  several  discrete  amplitude  levels  in  accordance  with  statistical  rules  to  be 
described.  Following  Skolnik,  *he  mathematical  treatment  will  include  not  only  the  important  case 
called  "Natural  Thinning",  in  which  the  aperture  element  density  varies  from  unity  (filled)  to  low 
(sparsely  filled)  in  an  attempt  to  approximate  the  desired  weighting,  but  will  also  consider  thinnirg 
beyond  this  "natural'  range.  The  density  of  elements  is  thus  scaled  by  a  constant  k  <  1  such  that  the 
probability  of  any  element  being  excited  can  be  <  k  (instead  of  unity)  while  the  ideal  aperture 
distribution  remains  proportional  to  the  naturally  thinned  case  (k  =  1).  Throughout  this  report,  the 
analysis  will  account  for  k  <  1.  but  the  examples  will  be  restricted  to  natural  thinning. 


5.  Taylor,  T.T.  (1960)  Design  of  Circular  Apertures  for  Narrow  Beamwidth  and  Low  Sidelobes.  IRE 
Trans.  AP-8:23-26. 
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A.  ARRAY  AMPLITUDE  TAPER  A(x,y) 


Figure  1.  Array  Amplitude  Taper 
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Figure  2.  (a)  Ideal  Taper  (Dashed)  and  Method  1  Source  Weight  Options  (b)  Ideal  Taper  (Dashed)  and  Method  2  Source  Weight 
Options  (c)  Ideal  Taper  (Dashed)  and  Method  3  Source  Weight  Options 
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2.1  Description  of  Methods 

2.]  1  METHOD  1 

Method  1  of  Figure  2a  shows  the  simplest  combination  of  thinning  with  amplitude  quantization: 
Discrete  weights  Vp  are  chosen,  always  above  the  Ideal  Illumination  envelope,  (except  at  the  origin) 
and  sources  are  either  left  "on"  or  turned  completely  "off1  in  accordance  with  an  algorithm  that  makes 
the  azimuthally  averaged  amplitude  at  any  radial  position  nearly  that  of  the  taper.  The  probability  of 
assigning  the  weight  Fn  =  Vp  to  an  element  at  location  pn  In  the  radial  annulus  pp_,  <  p  <  pp  is  given  by: 

P(Fn  =  Vp)  =  kAn/Vp  (1) 

where  A„  is  the  amplitude  of  the  ideal  illumination  at  the  nth  element.  The  probability  of  turning  the 
element  off  is  (k-P(Vp)].  The  figure  indicates  that  some  elements  are  left  on  at  the  value  Vp  while 
others  are  set  to  zero.  In  terms  of  a  distribution  around  some  annular  ring  of  elements,  the  above 
result  also  indicates  that  the  average  of  the  weights  around  a  large  ring  of  elements  will  also  be  given 
by  kA^/Vp.  Ii.  principle,  there  can  be  many  discrete  levels,  although  In  practice  some  small  number 
will  be  chosen. 


2.1.2  METHOD  2 

In  Method  2.  the  staircase  of  discrete  levels  is  again  kept  above,  or  Just  touching  the  ideal 
envelope  A^  but  in  this  case  the  average  level  for  a  ring  of  sources  is  attained  by  also  exciting  some 
fraction  of  elements  with  the  next  lower  level  of  amplitude.  Since  the  next  discrete  level  is  also 
available  to  a  solid  state  module,  it  would  seem  reasonable  to  take  advantage  of  this  added  degree  of 
freedom  instead  of  restricting  oneself  to  the  on-off  devices  In  Method  1.  For  Method  2,  whether 
considered  from  an  average  over  an  ensemble  of  arrays,  or  counted  from  the  annulus  of  a  large  array, 
if  there  are  M2  elements  at  the  level  Vp.  Mj  elements  at  the  next  lower  level  Vi  and  still  others 
missing  or  turned  off.  then  for  the  average  level  to  be  A„.  one  requires: 

M  2  M  j 

M  ,  +  M2  V?  +  M  |  +  M2  VP  +  i  =  An  • 


Since  the  total  excitation  probability  is  k,  to  account  for  generalized  thinning, 

P(Fn  =  Vp)  +  P(Fn=Vp+1)  =  k. 

We  associate  the  probability  of  exciting  elements  at  levels  Vp  and  Vp+1,  respectively,  with 


M. 


M, 


P(Fn  =  V  =  kM7Tl^’  P(Fn  =  Vp+1)  =  k 


(3) 

(4) 


Solving  for  the  probabilities  yields 

-  V  il  (Vp  -  An) 

P(F"  *  V  k  Vp  -  Vp, ,  •  P<Fn  =  V,)  =  k  v^7 


(5) 


The  sketch  illustrating  Method  3  looks  quite  different  from  that  of  Method  2,  but  in  fact  the  only 
geometrical  difference  is  in  the  selection  of  radii  pp.  In  Method  3.  the  staircase  of  discrete  levels  is  not 
kept  above  the  envelope,  but  is  roughly  half  below  and  half  above  the  taper  at  A(pp).  If  a  particular  pp  is 


i 
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selected*  the  step  height  is  chosen  so  that  it  is  bisected  by  the  A  (pp)*  that  is*  A  (pp)  =1/2  (Vp  +  Vp+i).  This 
choice  is  arbitraiy  and  can  lead  to  inferior  results  for  cases  with  very  few  steps*  as  will  be  discussed  in 
Section  4.4.  Method  3  must  contend  with  current  sources  in  circles  or  annuli  for  which  lies  either 
above  the  bisected  staircase  or  below  it.  The  consequence  is  that  two  sets  of  probability  densities  arc 
needed  to  implement  the  method*  as  will  now  be  discussed. 


2.1.3  METHOD  3 


In  Figure  3,  a  section  of  the  Method  3  configuration  has  been  expanded.  Sources  at  grid  locations 
denoted  by  ^  and  n2  are  displayed.  The  value  of  the  weight  at  nb  n2,  namely*  Fn  2  are  allowed  to  t 
either  of  two  values: 

Fm  =vp  or  Vp  +  1  :  Fn  2  =  Vp.  j  or  Vp  (6) 


The  probabilities  for  Fn  assuming  one  of  the  two  values  are  defined  below  by  the  "relative"  closeness  of 
A,,  to  a  level,  (See  Figure  3) 

Region  (I)  Step  above  taper 


|Anl  -  V  ,) 

P(Fn  ,  =  Vp)  =  k  ,v  -V  ^T  :  P(Fnl  =  W  =  k  -  P<Fn  1  =  Vp) 


(7) 


Region  (II) 


P(Fn2  =  Vp) 


IV,  “  An2l 
IV,  -  V  '' 


Step  below  taper 

P(Fn2=Vp.1)=k-P(Fn2  =  Vp) 


(8) 


Figure  3.  Expanded  View  of  Method  3  Geometry 
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Intuitively,  as  approaches  Vp>  the  probability  should  increase  that  the  weight  Fn  at  the  source 
location  pn  assumes  the  value  Vp.  Thus,  drawing  a  uniformly  distributed  random  variate  (rv)  from  a 
pseudo-random-number  generator  producing  samples  over  0  <  rv  <  k  lets  one  decide  on  the  value  given 
to  Fn.  For  example,  suppose  a  source  element  is  located  at  pn2*  a  pi  sition  where  the  taper  amplitude 
exceeds  the  level  value  (in  this  case  Vp).  If.  say  P(Fn2  =  Vp )  =  0.85  one  would  set  Fn2  =  VpM  for  rv  1 0.85. 
Using  the  probability  formulas  for  Methods  1  and  2,  respectively,  cue  can  decide  in  a  similar  fashion 
whether  to  set  Fn  =  Vp  or  0  and  Fn  =  Vp  or  Vp  +  1 .  The  thinning  algorithms  for  all  three  methods  were 
constructed  using  this  design  philosophy. 

2.2  Average  Power  Patterns  of  Thinned,  Quantized  Arrays 

If  the  ideal  aperture  taper  of  Figure  1  were  sampled  periodically,  and  radiated  from  a  filled, 
periodic  array,  the  sampled  values  An  would  radiate  to  produce  the  far  field  pattern 

N 

Eo(0.4>)=  Z  Anexp(jd>n)  (9) 

n=  1 

where  the  <I>n(0,<>)  is  the  phase  factor  as  detailed  in  Section  4.  In  the  case  of  the  new  array  with  a  finite 
set  of  amplitudes  and  the  source  configuration  subjected  to  statistical  thinning,  the  field  intensity 
pattern  is  written: 

N 

E(0,< 4>)=  z  Fnexp(JOn)  (10) 

n=  1 

where  the  Fn  are  selected  randomly  in  accord  with  the  previously  described  rules. 

Ensemble  averaging  the  field  intensity,  E(0,4>),  from  many  thinned  arrays  leads  one  to  set 

Fn  =  kAn  (11) 

in  order  to  have 

E(0,<t>)  =  kE0  (9.4))  (12) 

For  large  arrays,  we  assume  that  the  azimuthal  sample  average  of  the  Fn  is  nearly  the  same  as  the 
ensemble  average;  namely  that 

<  F  n  >  =  Fn  =  k A  n  .  (13) 

For  the  purpose  of  deriving  averaged  expressions,  the  ensemble  average  is  more  convenient  to  invoke. 
The  mean  value  of  the  field  strength  is  thus  given  by 

_  N  N  _ 

E(0.0)  =  I  Fnexp(jOn)  =  I  Fn  exp(j <hn)  (14) 

n=l  n=I 

=  kEo(0,<t>) 

which  is  the  pattern  of  the  ideal,  filled  array. 
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The  power  pattern  Is  the  absolute  square  of  the  field  Intensity  pattern.  Following  Skolnik,  the 
power  pattern  is: 

I  E(6,<J>)1 2  =  E(0.4>)  E»(e.4>)  (15) 


=  FmFnexP  U  ^ n )  1 


and 

I  E(0,4>)l 2  =  ZFn2  +  Z  SFmFnexp  [j(d>m  -  0>n)] 

m  *  n 

is  its  expected  value  for  real  Fn.  The  m  =  n  index  is  separated  out  to  invoke  the  theorem  from  statistics 
which  states  that  for  statistically  independent  random  variables,  the  mean  of  a  product  is  the  product 
of  the  means. 


2.2.1  METHOD  1 


At  this  point  the  analysis  departs  from  that  of  Skolnik  to  address  the  specific  algorithms  under 
study.  First,  it  is  convenient  to  break  up  the  sums  over  n  into  groups  that  correspond  to  the  discrete 
levels.  Thus,  a  sum  over  n  for  the  entire  array  is  written  as  a  double  sum  over  the  elements  at  all  P 
levels,  with  the  index  n  vaiying  between  Np_j  and  Np  in  the  pth  level  (that  level  between  Pp„1  and  Pp); 


£ 

n(p) 


(16) 


The  Quantities  Fn  can  assume  either  of  two  values,  Vp  or  O,  for  n  in  the  region  Np_i  <  n  <  Np.  Using 
the  definition  for  the  mth  moment  of  a  discrete  probability  density  function  (PDF), 

X™  =  Z  xm  P(x,  =  x],  (17) 


one  gets  from  Eq.  (1) 

Fj  =  Vp2(kAn/Vp)  +  0 


(18) 


for  i  =  1  and  2,  leading  to 

Z  F2  =  k  IVD  ZAn  .  (19) 

n  P  P  n(p) 

The  second  term  of  Eq.  (15)  is  evaluated  below  using  the  theorem  concerning  the  mean  of  the  product  of 
statistically  independent  variables. 

I  PmFnexpIj(Om-On)l  (20) 

m*n 

=  £  Z(FmFnexp[j(d>m-On)] 
rii  *n 


Using  Eq.  (13),  this  becomes 
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(21) 


I  E  k2AmAn  exp  U  COm  -  «I>n)) 
m*n 

=  k2 1  E0  (0,<J))l2  -  k2Z  (A n) 2 

n  11 

Combining  both  sums  to  produce  the  average  power  result  leads  to: 

I  E(e,<t>)  I  2=  k2 1 Eo(0,<t>)| 2  +  EV.E  kAn  (1  -  kAn/V_) 

P  v  n  (p)  ^ 

=  k2  |E0(9.4))I2  +  PSL  (22) 

where 

PSL=  EVpI  kAn(l  -  kAn/V  )  . 

P  *n  (p)  K 

The  expression  is  similar  to  Skolnik's  result  in  that  it  has  two  terms:  the  first  is  the  power  pattern  of 
the  ideal  filled  array,  and  the  second  is  a  constant  or  uniform  level  independent  of  the  angle  variables 
0  and  <t>.  It  is  a  generalization  of  Skolnik's  result,  and  reduces  exactly  to  it  if  there  is  only  a  single  level 
V0  and  not  a  set  of  discrete  levels.  Quantitative  results  from  the  method  will  be  discussed  in  a  later 
section. 


2.2.2  METHOD  2 

o 

Using  the  probabilities  given  in  Eq.  (5),  we  have  from  the  moment  formula  for  Fn: 

Tp  )2  -  (V  )2  VP+*1  +  ry  )2  -IVP  ~  Alll 

irnl  ~  lvpi  V  -  V  ,  +  'vp+l'  V  -  V  , 

r  vp  P+1  v  "  v"j-i 


VP  VP+1 


(23) 


Summing  over  all  n,  and  using  the  double  sum  notation  described  for  method  1,  we  arrive  at 

Nn 


Z(Fn)2  =  ^2  NZ  [An  (Vp  +  Vp  +  1)  -  Vp  Vp  +  1]. 


(24) 


"p-1 


The  second  term  of  the  average  power  pattern  is  again  given  by  Eq.  (21),  as  can  also  be  verified  by 
inserting  the  expression  for  Fn  into  this  term  and  using  the  probability  formulas  (Eq.  5)  for  Method  2. 
Collecting  terms,  we  get  the  final  result: 


I  E(0.q>)  I  2  =  k2!  Eo(0,()))l  2  -  Ek2(An)2 


Nn 


+k  upyA",v-tvP.l>-vPv1]. 

Clearly  this  is  in  the  form  of  Eq.  22,  with  the  average  sidelobe  level 


(25) 


Psh=  p  nfp)  (kAn  (VP  +  VP+l)  '  VP  VP  +  l1  "  n  k2(A")2 


(26) 
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2.2.3  METHOD  3 


The  probabilities  for  Method  3  have  been  defined  in  Eq.  (7)  and  (8).  Again,  the  ensemble  averaged 
power  pattern  is  given  by  E(0,  0)2  shown  in  Eq.  (15).  The  major  difference  between  Method  3  averaging 
and  the  others  is  that  now  one  has  to  split  the  sum  over  n  into  two  parts,  namely, 


£(Fn) 


s>  12  ^ 


,?(Fnl)2  + 


(1) 


.i)(F-)2 


(27) 


where  the  (1)  and  (2)  denote  the  two  regions  and  options  encountered  in  a  typical  interval  (pp_,  ,  pp). 
Using  Eqs.  (6)  and  (7)  for  the  probability  densities,  and  Eq.  (17)  for  the  first  moment,  we  obtain  for 
Region  1 


—  v  k(Anl  -  Vp+i)  k ( Vp  ~  A n  1 ) 

nl  ?  VP-VP  +  1  P+1  VP"VP  +  1 


For  Region  2, 


r~  v  k(V.-i-A„a)  k(An2-  V.) 

rn2  vp  V  i  -  V 
P-1  p 


'P-1  VP-1  -  vp 


(28) 


(29) 


Next,  we  will  need  to  consider  the  second  moments.  For  Region  (1). 


(F„i)2  =  (Vp)2P(Fn  ,  =  Vp)  +  (Vp  +  1)2P  (Fnl  =  Vp  +  1)  ; 
similarly,  for  Region  (2): 

(Fn2)2  =  (Vp)2P  (Fn 2  =  Vp)  +  (Vp. j)2 P(Fn2  =  Vp.!)  . 

Collecting  contributions  from  all  P  intervals,  namely,  for  the  entire  array,  gives 


(30) 


(31) 


I  E(0,p)l  =I(Fn)2+  T,  S  Fm  Fn  exp  j(<I>n  -  4>m) 


n  m 
m  *n 


(32) 


which  again  may  be  written  as 

=  Z  (Fn)2  +  k2  |eo(0.<J.)  |2  _  k2  E(An)2  (33) 

Substituting  the  probability  formulas  into  the  equations  simplifies  the  results  to 

I  E(0,p)l 2  =  k2 1  Eo(0 ,4> )l 2  -  k2  Z(A_)2 

n 

+  kp(mIAnl  (vp  +  vp-i)  -  vp  vp+i) 

+  (2)IA„2  (Vp  +  Vp-i)  "  Vp  Vp_,]J  (34) 

The  average  sidelobe  level  for  Method  3  is  thus  given  by 
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PSL  =  -k2  54!  ♦  kX  {l[A„,  <Vp  +  Vp„)  -  Vp  Vp.,| 


+  .|.|A„2  (Vp  +  vp.,)  -  Vp  Vp.,]j 


(35) 


3.  ARRAY  PARAMETERS 

The  several  expressions  obtained  above  give  average  patterns  that  have  two  terms;  the  Ideal 
filled  array  pattern,  and  an  average  pattern  Independent  of  angular  variation.  At  first  glance  it 
appears  from  these  expressions  that  one  need  only  evaluate  the  average  (far  sidelobe)  patterns  to 
compare  the  results  of  the  three  thinning  methods;  however,  there  are  other  array  parameters  that 
depend  upon  the  quantization  and  thinning.  They  are  evaluated  In  this  section  and  described  in  a 
number  of  figures  that  follow.  The  chosen  array  is  laid  out  on  a  rectangular  grid,  with  inter-element 
spaclngs  of  0.5  X  In  each  direction,  the  array  radius  for  the  pth  annular  ring  being  pp.  Within  each  ring 
there  are  a  fixed  number  of  elements,  and  the  several  methods  differ  In  the  way  they  perform  the 
statistical  thinning  and  quantizing  for  each  ring. 


3.1  Normalized  Input  Power 

If  there  are  a  total  of  N  elements  In  the  array  with  amplitudes  Fn,  then  the  total  power  into  the 
array,  normalized  to  the  Input  power  into  a  uniformly  Illuminated  filled  array  (Fn  =  1)  Is: 

E  (Fn)2 

P.n  =  — 53“—  •  (36) 


Pm  is  shown  in  Figures  9a,  10a,  1  la  for  the  three  methods. 


3.2  Average  Sidelobe  Levels 


The  (ensemble)  average  sidelobe  power  levels  PSl  given  in  Eqs.  (22),  (28)  and  (37)  should  be 
normalized  to  the  ensemble  average  peak  of  the  radiated  signal.  For  convenience  however,  it  has  been 
normalized  to  the  peak  of  the  specific  designed  array  or  (ZFn)2.  The  average  sidelobe  level  normalized 
to  the  beam  peak  is  given  below: 


SL  = 


PSL 

(£Fn)2 


(37) 


3.3  Array  Directivity 

The  directivity  of  an  array  depends  upon  the  element  pattern,  and  does  not  admit  to  a  closed 
form  expression  in  that  it  needs  to  Include  mutual  coupling  terms  explicitly;  otherwise  one  must 
resort  to  an  integration  of  the  entire  pattern.  Skolnik  uses  a  formula  that  neglects  mutual  coupling 
(and  assumes  a  unity  gain  element  pattern)  to  obtain: 
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(38) 


D  = 


2Fr 


2  F„ 


2 


2 


One  Interesting  point  about  Eq.  (38)  Is  that  it  is  not  sensitive  to  the  (ensemble)  average  sidelobe 
level  or  the  relative  degree  of  thinning.  For  a  single  discrete  level  with  each  Fn  either  one  or  zero,  the 
directivity  is  just  the  number  of  non-zero  biased  elements,  independent  of  the  array  size  or  degree  of 
thinning,  even  though  the  average  sidelobe  level  may  have  Increased  to  a  degree  where  it  should 
dominate  the  directivity.  Put  another  way,  this  formula  says  that  an  array  of  N  elements,  excited  with 
equal  amplitudes,  has  the  same  gain  for  any  array  spacing,  as  long  as  the  mutual  coupling  can  be 
neglected.  As  written,  the  formula  implies  a  unity  gain  (omnidirectional)  element  pattern,  and  so  o  ic 
can  argue  that  for  a  very  sparse  array  the  gain  should  be  multiplied  by  2  to  account  for  hemisphencal 
radiation  (while  still  neglecting  coupling).  When  the  array  elements  are  closely  spaced  (say  about 
0.5  X)  then  clearly  the  array  directivity  for  the  uniformly  excited  filled  array  has  to  approach  the  area 
directivity  (4rcA/X2),  and  one  can  show  this  as  being  consistent  with  Eq.  (38)  only  if  an  element 
directivity  of  n  is  assumed  to  multiply  Eq.  (38). 

One  can,  in  fact,  show  that  if  the  element  patterns  of  an  array  are  identical  (no  end  effects)  then 
the  array  directivity  can  be  written 


2F_ 


D  =  Dc 


2  F„ 


(39) 


and  the  only  problem  with  the  above  formula  is  that  the  element  pattern  De  varies  between  tc  for  a 
filled  array  lattice  (with  unexcited  elements  terminated)  and  2  for  an  array  with  hemispherical 
element  patterns. 

In  the  present  study,  we  wish  to  address  several  cases  in  which  the  array  is  nearly  filled,  and 
have  introduced  an  expression  for  the  average  directivity  of  a  statistically  thinned  array  radiating 
into  a  half-space.  For  all  three  methods  the  ensemble  average  pattern  of  the  array  takes  the  form: 

E(0.0)2  =  k2  |E0(e.4>)|2  +  PSL  (40) 


where  the  IEqIG,  4>)l 2  is  the  ideal  pattern  of  the  tapered,  filled  array.  The  directivity  of  an  antenna 
radiating  into  a  half- space  is  given  by 

4 n  Pm^y 

D  =  -  (41) 

P(0,4>)  d£2 


where  is  the  power  density  k2|  E^max  at  the  beam  peak,  £2  is  the  solid  angle  subtended  in  the  upper 
half-space,  and  d  £2  =  sin  0  d  0d  q>. 

The  average  directivity  in  this  case  is: 


J 


4tc 


k  Eo(0,4>)  2dQ  + 


si. 


J 


d£2 


1  +  1/2  k-2  D0  SL 


(42) 


for  the  sidelobe  level  SL  =  PS(/Pinax  and  D0  is  the  directivity  of  the  ideal  pattern. 
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Eq.  (42)  shows  that  the  average  directivity  depends  on  the  ideal  pattern  directivity  and  the 
average  sidelobe  level.  When  the  average  sidelobe  level  approaches  zero,  as  for  the  filled  ideal  aperture 
taper,  the  directivity  tends  to  the  correct  limit  [of  4n(AA2eap]  where  eap  is  the  aperture  efficiency  of  the 
chosen  distribution  (in  this  case  the  circular  aperture  Taylor  distribution).  When  the  average 
sidelobes  are  at  the  isotropic  level,  the  directivity  is  reduced  by  about  1.8  dB  (SL  -  1/Dq).  For  very  high 
average  sidelobes,  the  directivity  approaches  2/SL  above  the  sidelobe  level,  which  is  the  case  when 
virtually  all  the  power  is  in  the  average  sidelobes  and  the  radiation  is  contained  in  a  half-space. 

Notice  that  the  thinning  constant  k  does  not  enter  explicitly  in  this  formula,  but  is  included  in  a 
sidelobe  level. 

The  expression  (42)  is  not  exact,  for  it  implies  an  omnidirectional  (half-space)  element  pattern 
applied  to  the  integration  of  the  average  sidelobe  power  in  the  dominator  of  Eq.  (41),  and  so  neglects 
element  pattern  changes  due  to  coupling.  If  it  is  assumed  that  the  array  elements  are  closely  spaced,  a 
cos  0  element  pattern  might  be  assumed  to  modify  the  average  sidelobe  power  integral.  That  would 
lead  to  a  factor  of  1  /4  multiplying  the  term  SL  for  this  case  of  very  little  thinning.  However,  in  that 
limit  the  average  sidelobes  are  very  low  and  the  directivity  very  nearly  the  ideal  (so  it  doesn't  matter 
what  element  pattern  is  used).  Alternatively,  when  the  array  is  highly  thinned,  the  element  pattern  is 
more  nearly  omnidirectional  (in  the  hemisphere),  and  expression  (42)  should  then  be  a  good 
approximation.  The  validity  of  it  was  tested  by  integrating  the  directivity  formula  (41)  assuming  q> 
symmetry,  for  several  cases  as  noted  on  the  figures. 

3.4  Normalized  EIRP  (Effective  Isotropic  Radiated  Power) 

The  product  of  normalized  input  power  and  array  (average)  gain  is  tin  Effective  Isotropic 
Radiated  Power  (or  EIRP).  The  EIRP  is  one  of  the  primary  measures  of  the  antenna  system,  as  it 
produces  the  power  density  at  the  target  (or  receiver).  For  a  lossless  antenna  system  the  directivity  and 
gain  are  equal  and, 

EIRP  =  D  Pln  (43) 

4.  ARRAY  CALCULATIONS  AND  RESULTS 
4.1  Array  Pattern  Calculations 

The  array  factor  for  a  circular  aperture  array  with  elements  located  on  a  rectangular  grid  is: 

E(0,4>)  =  K I E  Fm  n  exp  {  J(2tc/X)  (md*  (u  -  Uq)  +  ndy  (v  -  Vq))  ) 
for 


u  =  sin  0  cos  <)>;  u0  =  sin  0O  cos 

v  =  sin  0  sin  <J>;  v0  =  sin  0O  sin  4>0  (44) 

for  m  and  n  within  the  circular  planar  region  as  shown  in  Figure  1,  and  elements  centered  at  the 
locations  (x^  yn)  =  (mdx.  ndy). 

The  constants  Fmn  are  real  and  are  determined  from  the  statistical  algorithms  given  in  previous 
sections.  The  constant  K  is  chosen  to  produce  max  F(0,  <J>)  =  1;  namely 
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Since  the  array  Investigated  in  this  paper  contained  several  thousand  elements,  the  circular  Taylor 
amplitude  distribution  [Ref  5[  (characterized  by  the  parameter  n )  had  to  be  sampled  thousands  of  time*  at 
the  nodes  of  a  rectangular  grid  locating  the  current  elements.  To  speed  the  process  of  evaluating  the 
Taylor  function,  we  sampled  it  at  19  positions  across  its  radius  and  used  a  cubic  spline  code  for  rapid 
interpolation  at  the  grid  radii,  r  (xm,  yn),  leading  to  significant  reductions  in  program  execution  time. 

The  Taylor  illumination  being  emulated  was  always  chosen  with  n  as  a  large  as  possible  while 
maintaining  monotonically  decreasing  current  distributions  from  the  origin.  This  selection  leads  to 
the  highest  aperture  efficiency  and  narrowest  beamwidth  subject  to  given  sidelobe  criteria  while 
maintaining  realistic  taper  illuminations.  The  selected  values  of  n  accompanied  by  their  compiutd 
aperture  efficiencies  are  given  below: 


Level  (dB) 

n 

Efficiency 

-30 

5 

0.8623 

-35 

6 

0.7880 

-40 

7 

0.7186 

-45 

9 

0.6620 

-50 

11 

0.6106 

-55 

13 

0.5651 

-60 

15 

0.5251 

-65 

18 

0.4909 

-70 

21 

0.4603 

-75 

23 

0.4325 

These  efficiencies  were  calculated  from  Hansen6  using 


F  2 

1  n 


nti  Jo2  (»Hn) 


(46) 


where  the  pn  are  the  zeros  of  J  i(rcp). 

These  efficiencies  were  used  in  the  computation  of  directivity  for  the  ideal  pattern  using: 
4xA 


D0  = 


A.2 


(47) 


The  pattern  calculations  for  a  filled,  sampled  aperture  of  this  size  are  very  nearly  identical  to 
those  of  the  continuous  aperture  with  the  Taylor  taper  and  so  exhibit  sidelobes  at  the  design  level. 


4.2  The  Design  Process 

Since  the  objective  is  to  apply  statistical  thinning  methods  to  a  sampled,  step-wise  continuous 
aperture  illumination,  the  following  question  arises:  What  amplitude  levels  (Vp)  and  radii  (pp)  are  to 
be  assigned  to  the  steps?  Ideally,  according  to  some  criteria,  the  best  (Vp,  pp)  is  determined  and  the 


6.  Hansen,  R.C.  (1983)  The  Handbook  of  Antenna  Design,  Vol  2,  Chapter  10,  pp.  157,  Rudge  et.  al. 
(Editors),  Peter  Peregrinius  Ltd. 


array  is  thinned,  simultaneously.  The  inherent  difficulties  of  global  optimization  in  this  context 
leads  us  to  the  (suboptlmal)  design  scheme  described  below.  In  it,  the  best  (pp)  is  determined  first,  and 
then  the  array  is  thinned. 

We  assume  that  by  approximating  the  ideal  taper  with  the  staircase  function  shown  in  Figure  lc. 
Taylor  power  pattern  features  will  emerge  as  the  number  of  steps  (P)  increases  indefinitely.  To  choose 
the  "best"  power  pattern  for  a  small  number  (P)  of  steps  we  proceed  as  follows: 

1.  N  sets  of  randomly  generated  pp  are  created; 

Corresponding  sets  of  Vp  are  obtained  by  evaluating  the  Taylor  taper  at  the  pp: 

3.  The  power  patterns  for  all  N  staircase  apoaizations  are  calculated: 

4.  The  interval  containing  the  first  few  Taylor  sidelobes  (in  each  power  pattern)  is  sampled  and 
the  number  of  ordinates  (Q)  exceeding  the  design  sidelobe  level  (SLL)  is  counted:  and  finally, 

5.  The  staircase  parameters  (pp,Vp)  yielding  the  smallest  Q  is  deemed  the  best  for  the  subsequent 
thinning  operations. 

Typically,  N  =  200  configurations  were  considered. 

The  power  patterns  required  In  Step  3  are  easily  calculated  using  the  two-dimensional  Fourier 
transform  over  the  circular  aperture  with  constant  illumination  in  each  ring.  By  virtue  of  the 
rotational  symmetry,  one  is  led  to  the  evaluation  of  a  finite  Hankel  transform  for  each  ring.  For  P 
steps  in  the  staircase,  the  far-field  amplitude  distribution  becomes 


Fs(sin0)  =  CI[Vp-Vp  +  1)(Pp)2 


( J  i(2xpp  sin  6)} 
2xpp  sin  9 


(48) 


with  C  a  normalization  and  pp=  rp/X,  Vp+1  =  0,  Fs2  gives  the  required  power  pattern. 

Some  of  the  best  patterns,  obtained  for  various  specifications  of  desired  sidelobe  level,  were 
those  for  which  the  radii  varied  as 


Pp  -  PAP  p=  1,2,3... 

for  Methods  1  and  2.  For  Method  3,  no  simple  relationship  governing  the  radii  could  be  inferred. 

Having  obtained  the  implied  pp,  the  three  algorithms  were  run  to  evaluate  the  patterns  as  well  as 
to  determine  how  the  basic  parameters  -  gain,  average  sidelobe  level  and  normalized  EIRP  -  varied 
with  method.  The  results  are  summarized  in  Figures  9  through  14, 


4.3  Array  Dimensions  and  Generalizations 

In  nearly  all  cases  the  data  was  taken  for  the  following  array  dimensions: 

Array  Diameter:  51  K 

Element  Spacing:  0.5  \  (square  grid) 

No.  of  Elements  in  Filled  Array:  7845 
Design  Sidelobe  Level  (dB):  -30  to  -75 

Although  most  data  was  collected  for  only  one  array,  it  is  possible  to  infer  the  average  parameters  of 
other  large  arrays.  It  is  shown  in  the  Appendix  that  for  a  large  array  with  a  given  taper,  the 
normalized  sidelobe  level  varies  as  1/N,  where  N  is  the  total  number  of  array  elements.  The  array 
directivity  varies  directly  with  N.  This  means  that  the  product  N  SL  and  the  quotient  D/N  are 
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constants,  independent  of  array  size,  once  the  arrays  are  sufficiently  large.  For  this  reason  the  curves 
in  Figures  9  through  14  are  generalized  with  ordinates  on  the  left  that  include  the  array  element  count 
N,  while  the  ordinates  at  the  right  of  each  figure  are  the  actual  computed  data.  The  ordinates  of  the 
sidelobe  curves  are  given  a  normalized  sidelobe  level  that  is  the  product  of  N  times  the  sidelobe  ratio. 

10  log jq  SLnorm  =  10  log10  SL  +  10  logi0N 

The  ordinate  of  the  normalized  E1RP  curve  is: 

10  log10  Dnorm  =  10  log10  D  -  10  log10N 

The  ordinate  of  the  normalized  EIRP  curve  is: 

10  log]0  (EIRPnonn)  =  10  logj0  EIRP^^  -  20  logj0N 

Since  N  was  7845  elements,  the  parameter  10  log  N  is  +  39  dB. 

4.4  Discussion  of  Results 

Figure  4a  shows  the  power  pattern  resulting  from  a  discrete  array  of  elements  located  at  the 
nodes  of  a  rectangular  sampling  grid  (x^y^  Ax  =  Ay  =  0.52)  within  the  circular  aperture.  The  grid 
samples  a  -50dB  Taylor  taper,  with  amplitude  A  (rmn)  for 

(rnvn  =  ^4  +  y„  )l 

the  pattern  is  clearly  a  good  approximation  to  that  arising  from  the  unsampled  A(r). 

Figure  4b  shows  the  pattern  of  a  filled  array  with  three  discrete  amplitude  levels  used  to 
approximate  a  -50dB  Taylor  illumination  (see  Figure  lc).  The  amplitude  levels  (w)  and  quantization 
step  radii  (p)  are  shown  in  this  Figure,  as  well  as  in  those  that  follow.  Since  the  array  is  not  thinned 
there  are  some  relatively  high  sidelobes,  but  the  annular  distribution  of  the  discrete  amplitude  levels 
tends  to  minimize  grating  lobes.  The  goal  of  this  study  is  to  improve  this  and  similar  patterns  by 
utilizing  various  statistical  thinning  procedures  in  addition  to  the  discrete  level  approximation  to  the 
taper. 

Figure  5a  shows  the  result  of  thinning  the  array  using  the  same  excitation  amplitude  at  all 
elements  and  statistically  thinning  the  array  (as  done  by  Skolnik)  to  match  the  -50  dB  Taylor  taper 
illumination.  Clearly  the  pattern  has  the  correct  beam  width,  but  has  many  sidelobes  above  the  -35dB 
level  and  average  sidelobes  at  -37.5  dB.  The  pattern  directivity  is  36.2  dB.  The  degree  of  thinning  is 
indicated  in  Figure  5b,  which  shows  a  continuous  grading  of  elements  removed  to  approximate  the 
desired  density  taper.  In  this  figure  the  dashed  lines  indicate  elements  that  are  not  excited.  Notice 
that  the  attempt  to  match  an  exceedingly  low  sidelobe  taper  has  led  to  very  inefficient  use  of  the 
aperture  and  to  high  average  sidelobes.  This  is  a  rather  extreme  example,  because  if  the  array  had  been 
larger  (by  a  factor  of  10  for  example)  then  the  average  sidelobes  in  Figure  5a  would  be  reduced  another 
10  dB.  As  shown  by  Skolnik,  the  average  sidelobe  level  is  approximately  1/N  for  highly  thinned 
arrays.  Figures  6,  7,  and  8  show  the  same  example  with  three  steps  of  amplitude  illumination.  For 
comparison,  recall  that  Figure  4b  shows  the  optimized  three-step-pattem  without  thinning  and 
indicates  peak  sidelobes  in  the  neighborhood  of  ~27  dB.  In  7b  and  8b  dashed  lines  indicate  that  the 
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(a)  Power  Pattern  from  -50  dB  Taylor  Taper  Sampled  for  dx  =  dy  =  \/2 


(b)  Power  Pattern  from  a  Filled  (Unthinned),  3-Step  Distribution 


Figure  4. 


(b)  Distribution  of  Non-Excited  (Thinned)  Elements  for  Figure  5a 

Figure  5. 
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(a)  Power  Pattern;  3-Step  Distribution;  Method  1 


(b)  Distribution  of  Non-Excited  (Thinned)  Elements  for  Figure  6a 


(a)  Power  Pattern;  3-Step  Distribution;  Method  2 


(a)  Power  Pattern;  3-Step  Distribution;  Method  3 


KIH3*  3.0  HOT  ICC  THAT  H  *54  15  A$$l*£0 


(b)  Distribution  of  Elements  with  Increased  (+)  or  Decreased  (-)  Excitation  for  Figure  8a 


element  is  excited  at  the  level  of  the  next  (lower)  step,  and  in  8b  the  plus  signs  indicate  that  the 
previous  (higher)  step  level  is  used. 

Figure  6a,  7a,  and  8a  show  the  results  of  applying  the  three  algorithms  to  smooth  the  average 
amplitude  illuminations.  Apparently  methods  2  and  3  have  very  similar  results,  with  average 
sidelobe  levels  between  -45  (Method  3)  and  -47  dB  (Method  2).  Method  1  has,  as  expected,  the  highest 
sidelobes,  at  -42.5  dB  but  these  are  substantially  suppressed  compared  to  the  case  of  a  single  step 
(thinning  with  uniform  amplitude)  (Figure  5a).  Peak  sidelobes  along  the  cut  <J>  =  0  range  between  -33  db 
(Method  1)  to  -38  dB  (Method  2). 

Figures  6b,  7b,  8b  show  the  state  of  the  array  aperture  for  the  three  methods.  Figures  9  through  1 1 
compare  the  results  from  the  3  synthesis  approaches  for  a  range  of  design  sidelobe  levels,  and  reveal 
some  interesting  differences  among  the  methods.  Figure  12  gives  the  directivity  and  selected  has  a 
function  of  design  sidelobe  level.  Figures  9,  10,  and  1 1  each  have  three  parts,  and  show  normalized 
input  power,  directivity  and  average  sidelobe  level  as  a  function  of  design  sidelobe  level  for  the  three 
methods.  Each  figure  shows  data  for  one,  two,  four  and  six  steps. 

For  one  step,  all  of  the  data  are  the  same;  the  three  methods  are  identical,  and  equivalent  to  the 
Skolnik  technique.  As  more  steps  are  added,  the  three  approaches  become  characteristically 
different.  In  almost  all  cases  the  input  power  and  directivity  decrease  with  decreased  design  sidelobes 
(as  one  would  expect  because  of  increased  thinning)  and  the  average  sidelobe  level  Increases  for  the 
same  reason. 

Figure  9  shows  the  directivity  of  an  array  synthesized  by  Method  1.  Several  additional  data 
points  have  been  added  to  this  figure  to  compare  the  several  directivity  formulas  Eq.  (38)  and  Eq.  (41) 
with  the  exact  gain  as  computed  by  numerically  performing  the  Integral  in  Eq.  (41).  Two  cases, 
representing  design  sidelobes  of  -30  dB  and  -60  dB  were  investigated  for  the  array  with  a  single  step 
height,  and  a  third  case  was  evaluated  for  design  sidelobes  of  -60  dB  and  six  steps.  The  exact  data  are 
given  by  the  circles,  which  are  shown  to  be  less  than  0.2  dB  from  the  plotted  results  using  the  formula 
(42).  The  results  of  using  Eq.  (39)  are  denoted  by  an  asterisk.  For  one  step  and  assuming  element 
directivity  of  n,  Eq.  (39)  leads  to  directivities  between  0.9  and  1.0  dB  higher  than  the  exact  number.  For 
the  case  with  one  step,  identical  accuracies  are  obtained  using  Methods  1  and  2.  For  an  array  with  six 
steps  and  Method  1  thinning,  the  results  are  closer,  and  here  the  results  of  using  the  formula  (41)  are 
within  0.3  dB  of  the  correct  value.  Figure  10,  computed  using  Method  2  thinning,  shows  only  the 
integrated  value  (at  -60  dB  design  sidelobe  level  and  six  steps).  In  this  case  all  three  methods  yield 
results  within  0.05  dB  of  each  other.  These  results  strongly  support  the  use  of  Eq.  (42),  which  is  used 
throughout  Figures  9,  10,  and  11. 

For  Method  1  (Thinning),  the  average  sidelobe  reduction  due  to  increasing  the  number  of  steps 
from  one  to  six  is  about  8  dB,  independent  of  the  design  sidelobe  level.  This  is  less  true  of  Methods  2 
and  3,  both  of  which  exhibit  a  wider  spread  of  average  sidelobe  level  as  a  function  of  the  number  of 
steps. 

There  is  an  interesting  anomaly  in  the  Method  3  data  for  two  steps.  As  shown  in  Figure  1 1,  the 
two-steps  curve  appears  quite  different  from  the  curves  for  1,  4,  and  6  steps.  The  reason  is  that  for 
Method  3  one  must  select  both  the  step  radius  pn  and  the  step  height  Vp  at  each  radius  to  define  the 
stepped  distribution.  The  method  chosen  for  determining  step  height  was  to  have  the  taper  amplitude 
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AVERAGE  SL  +  10  loglrt  N  DIRECTIVITY  -10  logirt  N  NORMALIZED 


Figure  9.  Pin,  Directivity,  and  Average  SL  vs  Design  SL;  Method  1 
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AVERAGE  SL  DIRECTIVITY  (dB) 


AVERAGE  SL.+ 10  log  N  DIRECTIVITY  -10  log  N  NORMALIZED  P 


DESIGN  SL„D 

0  D 


Figure  10.  Pin,  Directivity,  and  Average  SL  vs  Design  SL;  Method  2 
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AVERAGE  SLjd  DIRECTIVITY  (dB) 


Figure  12.  Taylor  Directivity  and  Selected  n  vs  Design  SL 

at  the  step  radius  equal  the  average  value  of  two  adjacent  levels,  namely, 
A(Pp)=^VP  +  vP  +  ll> 


Vp  -  A(Pp)  =  A(Pp)  -  Vp+  j  .  (49) 

Selection  of  this  condition  was  quite  arbitrary,  and  works  when  there  are  enough  steps.  However, 
when  there  were  only  two  steps  and  with  low  sidelobe  aperture  illuminations,  this  choice  resulted  in 
selection  of  a  second  step  V2  that  would  have  gone  below  zero  for  design  sidelobe  levels  less  than 
-50  dB.  To  avoid  this  outcome,  the  height  V2  was  set  to  a  very  small  positive  number.  The  net  result, 
shown  in  Figure  1 1.  is  for  the  two  step  data  of  Method  3  to  become  essentially  equal  to  the  one-step 
data  when  the  sidelobes  become  lower  than  -50  dB.  This  result  occurs  only  for  the  two-step  case,  and 
indicates  that  one  must  exercise  some  care  in  choosing  the  step-heights  for  Method  3.  Figure  12  shows 
the  normalized  directivity  of  the  ideal  Taylor  pattern  and  is  to  be  compared  with  Figures  9,  10  and  11. 
Note  that  several  of  the  sLx-step  configurations  for  Methods  2  and  3  have  nearly  the  maximum 
available  directivity.  Figure  12  also  shows  the  selected  n. 

Figures  13a  and  13b  give  additional  perspective  on  array  design.  The  figures  show  the  variation 
of  normalized  EIRP  with  design  sidelobe  level  for  Methods  1  and  2  with  one  to  six  steps.  Method  3 
results,  though  not  shown,  are  virtually  the  same  as  for  Method  2.  The  remarkable  fact  illustrated  is 
the  near  equality  of  EIRP  for  all  methods  at  any  given  sidelobe  level.  The  net  spread  of  these  curves  is 
approximately  1  dB  in  all  cases.  So,  in  terms  of  EIRP,  only  the  design  sidelobe  level  is  important,  not 
the  method  used  or  the  number  of  quantization  levels  (especially  for  greater  than  one  level). 
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Figure  13.  Normalized  EIRP  vs  Design  SL  for  Methods  1,  2  with  1-6  Steps 
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DIRECTIVITY  P  *  DIRECTIVITY 


Figures  9c  through  1  lc  show  average  sidelobe  levels  as  a  function  of  the  design  sldelobe  level  for 
various  numbers  of  steps.  The  plotted  sidelobe  level  parameter  is  [Avg  si]  dB  +  10  log10  N,  and  so  the 
data  is  invariant  with  the  number  of  filled  array  elements  (once  the  array  exceeds  a  minimum  size) 
However,  it  is  meaningless  to  select  an  illumination  to  produce  (near)  design  sidelobes  lower  than  the 
average  sidelobe  level,  and  so  it  is  of  interest  to  determine  how  large  an  array  must  be  in  order  to 
sustain  any  given  near  and  far  sidelobe  level.  A  good  estimate  of  this  quantity  is  obtained  by  setting 
the  average  sidelobe  level  equal  to  the  design  (peak)  sideiobe  level  minus  some  margin.  The  offset,  or 
margin,  accounts  for  the  fact  that  the  peak  sidelobes  of  the  random  distribution  may  be  significantly 
above  the  statistical  average  level,  and  it  is  common  practice  to  assume  a  margin  of  10  to  20  dB 
depending  on  the  percentage  of  sidelobes  allowed  to  exceed  the  peak  sidelobe  estimate. 

Using  the  normalized  sidelobe  data  from  Figures  9  and  10,  Figures  14a  and  14b  are  obtained  by 
setting 

SLjb  =  SLnorm  -  10  logjo  N 


-  ®  ^design  “  Margin^  (50) 

for  Margin  dB  a  positive  number.  Solving  for  N  yields: 

10  log10  N  =  SLnorm  -  SLdcs)gn  +  Margin^  (51) 

To  give  the  required  array  size,  the  figures  are  plotted  for  Methods  1  and  2  with  design  sidelobes  (the 
peak  Taylor  sidelobes)  set  to  values  between  -30  and  -75  dB,  and  with  Margin  dB  set  to  zero.  Thus, 
Figure  14  gives  the  minimum  array  size  required  for  the  average  sidelobe  level  to  equal  the  design 
near-sidelobe  level.  Note  that  for  P  =  1  (one  level)  these  cases  reduce  to  the  Skolnlk  et  al  result,  and  as 
shown  in  the  reference,  the  average  sidelobe  level  should  be  approximately  equal  to  the  number  N  of 
elements  remaining  in  the  highly  thinned  array.  Figure  14  confirms  the  results  for  the  arrays  studied 
herein.  The  figure  can  be  used  for  a  non-zero  margin  by  simply  increasing  N  by  the  value  of  the  desired 
margin,  for  example,  by  the  factor  10  for  each  desired  10  dB  of  margin. 

5.  CONCLUSION 

This  report  has  investigated  the  use  of  statistical  thinning  and  quantized  element  weights  to 
produce  low  sidelobe  patterns  using  large  circular  array  apertures.  The  major  results  of  the  analysis 
show  that  by  using  either  of  several  thinning  algorithms  it  is  possible  to  obtain  substantial  sidelobe 
reduction  by  a  combination  of  statistical  thinning  and  the  use  of  discrete  amplitude  quantization. 
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-35  -  40  -  45  -  50  -  55  -  60  -  65  -  70 

DESIGN  SIDELOBE  LEVEL  (dB) 

ARRAY  ELEMENT  COUNT  N 

(a)  Number  of  Array  Elements  for  Average  SL  =  Peak  SL;  Method  1 


ARRAY  ELEMENT  COUNT  N 

(b)  Number  of  Array  Elements  for  Average  SL  =  Peak  SL;  Method  2 


Figure  14. 
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Appendix 

Variation  of  Average  Sidelobe  Level  with  Element  Count  (k  =  1  Case) 


It  is  of  interest  to  scale  the  diameter  of  the  aperture  and,  correspondingly,  the  radial  dependence 
of  the  staircase  function  and  the  taper,  while  retaining  the  same  X/2  model  mesh  spacing.  To  effect 
this  scaling,  let 

D  — >  cD , 


rp->  crp. 


A(r)  A(cr), 
dx  =  dy  =  1/2 


(Al) 


(X/2  spacing). 


and 


{Vp  :  p  =  1 ,  2 . P}  remains  unchanged. 

As  noted  previously,  the  normalized  average  sidelobe  level  for  Method  2  is  given  by 


SL  = 


Pc 


I  I  (An(Vp+Vp  +  1)-VpVp+1)  -  £A„ 
i=l  n(p)  t  J  n 


SL 


(P-) 


(A2) 
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The  denominator  is  random,  but  even  for  "moderately  large"  arrays  its  ensemble  average  is  nearly 
that  from  a  single  spatial  realization  resulting  from  Method  2  thinning.  Thus, 


(  !Fn)2  = 

( LFn)2 

In  / 

V  n  / 

any  sample 
array  realization 


The  first  term  in  the  numerator  of  SL  is  given  by 


An(Vp+Vp+1)-V. 


pVi 


}• 


(A3) 


(A  4) 


Taking  account  of  the  radial  dependance  of  A^  and  interpreting  S[  dxdy  as  an  area  calculation  using  a 
two-dimensional  rectangular  quadrature,  then  transforming  to  polar  coordinates  yields 

2k  a 

Sj=(dxdy)_1  J  d<J>  Jdr  r  |A(r)  v(r)  -  v(r)  j 
oo 

1 

=  (dT^r)  a2  J  dpp  {A(p)v(p)-'0(p)}  (a 5) 

0 

P2 

=  N'2  ■  2k  =  {(vp  +  vp  +  1)  J  dp  p  ACp)  -  Vp  ^Pil  jpp  -  Pp-1  j  j 

P_  Pp-l 


=  N'% 

where  N'  =  p0  -  0  and  pi  =  1.  The  integral  of  pA  (p)  can  be  evaluated  in  closed  form  as  follows. 
Substituting 


n-  1 

A  ( p)  =  £  DmJ0(7Cgmp) 

m=0 

into  the  integral  yields 


(A  6) 


f  n  “  1  D  r  1 

dp  p  A(p)  =  I  -f“  jPpJj  UHmPp)- Pp-l  Jl  Pp-l)[- 
J  m=0  y  ' 


*‘P-  1 

This  completes  the  determination  of  ?.p  and  the  term  Sj. 
As  shown  above,  the  numerator  of  SL  can  be  expressed  as 

>-pN'2-I  A*, 
n 


(A  7) 
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where 


^P=  (pPl.  n,  P,  ...) 


but  not  N'.  We  are  left  with  the  evaluation  of 


n 


(p»r 


For  a  circular  Taylor  taper. 


1  f  r  \ 

A(r)=m?0  °m  J°  i7t^ma  J 


where  0  <,  r  <  a  =  N'  dx .  Letting,  P  =  r/a  and  approximating  the  area  under  A2(r)  with  a  rectangular 
quadrature,  one  obtains 


J  d(}>  Jdr  r  A2(r)  =  dx  dy  £  A^ 


(A  10) 


which  reduces  to 


v  a2  2* 
J  An  =  dx  d 


— y 
dx  dy  m= 


lo  J" 


dp  p  JQ  (TC^p) 


by  virtue  of  orthogonality.  The  integral  is  readily  evaluated  in  closed  form.  (See  Van  Bladel.  p.  514, 
(23)).  leading  to 

X  An  =  nN  £  Dm  Jo 
n  m=0 


=  N'2  C2. 
Thus  far,  one  obtains 

*-PN'2  -  CoN'2 


(All) 


(A  12) 


(  T .  c  1  2 
V  “*n; 
n 


To  evaluate  the  denominator,  we  note  that 


35 


(p") 


=  £  Fn  Fn' 


n,  n 


:  ,  Fn  Fn' 

n,  n 


Z  F  F  '  +  T  F  2 

^  1  n  1  n  ^  /L  1  \\ 

n 


n,  n 
n*n 


=  Fn  Fn'  +  Z  Fn2 

n,  n  n 


(A  13) 


=  £  AnAn'  +  I  2  {a„(v p+V p+1)-V pV p+1} 

n.n  P  n(p)  1  K  1 


n*n 


=  £  AnAn'+XPN 


,2 


n,n 

n*r/ 


=  Z  An  An'  -  X  An  +  A.PN 
all,  n 

n,  n 


=  numerator  of  SL 


It  follows  that 
SL  =  — 


1  + 


£  An  An' 


all 

n,n' 


X.pN'2  -  C2N'2 


(A  14) 


Next,  we  must  look  closely  at  the  sum 
S2=£  AnAn', 


(A  15) 


all, 
n.  n 


and  it  will  be  useful  to  revert  back  to  the  original  2-dimensional  array  notation,  namely. 


A  a  >  A„  ,,, 


A n’  ~>  A  n' m'. 

The  sum  may  now  be  written  as 


(A  16) 
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(A  17) 


$2-^  An  rn  An,ni 

all  , 
n,  n  , 
ni  m 


- S  X  X  S  An,m  An  +  (i  m  +  )  , 
n  m  i  i' 

which  is  a  2-dimensional  autocorrelation  function. 

Because  of  the  circular  symmetry  of  the  Taylor  function  and  the  circular  aperture,  it  is  useful  to 
convert  the  expansion  to  polar  form.  It  is  clear  that  the  2-dimensional  displacement  variables,  i  and 
i'  must  also  be  transformed  to  polar  form.  The  "area"  is.  therefore,  given  by 


2k 


2k 


2  2 

dx  ^v  ^2  ~ 


J  d9  J  dxx  J  d<()  J  dr  r 

o  o 


A(r)  A{r  +  x) 


The  change  to  dimensionless  variables 
p  =  r/a  and  u  =  x/a 
yields 


(A  18) 


(dx  dy)2  S2  *  (2ti)2  (N'  dx)4 


J  do  u  J  dp  p  A (p)  A(p+o) 

o  o 


=  (2jt)2  (N'  dx)4  I. 

Since  dx  =  dy  has  been  assumed,  and  a  =  N'  dx.  S2  may  be  written  as: 
S2  =  4it2 1  N'4. 

Inserting  into  the  expression  for  SL,  we  are  led  to 
1  1 


SL  = 


1  + 


IS  - 


1  +  Cs  N 


71 


(A19) 


(A20) 


(A21) 


(irfe)" 


.  since  N  =  JtN 


,2 


1  + 


Numerical  studies  show  that  to  an  excellent  approximation 
SL~  O  . 


r  a  ooi 


where  N  is  total  number  of  nodes  within  (and  on)  the  circular  aperture. 

The  evaluation  of  the  integral,  I,  can  be  accomplished  efficiently  and  accurately  by  numerical 
means.  Both  the  inner  and  outer  integrals  lend  themselves  well  to  low-order  Gauss  quadrature  cf 
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moments,.  Furthermore,  A(p)  can  be  rapidly  evaluated  using  the  same  cubic  snllne  code  that  was 
discussed  earlier  for  the  "thinning"  code.  It  should  be  apparent  that  Cs  (I,  Xp,  C2)  needs  to  be  determined 
only  when  p,  *  iAp  is  of  interest.  The  results  can  be  applied  to  arbitrary  aperture  diameters  as  long  as 
p,=  iAp. 
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